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ABSTRACT 

Approximate  methods  for  predicting  *e  center  of 
pressure  of  multi-segment  wmgs  are  presente  .  J 
5J  with  breaks  in  the  leading  edge  sweep  are 
considered.  The  methods  transform  a  multi-segjnent 
planform  into  an  equivalent  straight  tapered  P^oim 
Comparisons  are  made  with  predictions  from  a  vortex 
lattice  code  and  experimental  data  at  subsonic  speeds  and 

sunersonic  speeds.  The  most  accurate  method  studied  had 
an  average  error  in  predicted  center  of  pressure  location 
of  approximately  two  percent  of  the  root  chord. 

T  1ST  OF  SYMBOLS 

cc  angle  of  attack 
AR  aspect  ratio,  b2/S 

Cm  Etching  moment  coef.,  (pitching  moment/q-Scr) 
CN  normal  force  coef.,  (normal  force/qJS) 
cr  root  chord 

c,  tip  chord 

A  sweep  angle  m 

q  lateral  position  of  break  m  sweep  angle,  y/(b/Z) 

X  taper  ratio  (cjct) 

qTO  dynamic  pressure 

RN  Reynold’s  number 

S  wing  area 

.  center  of  pressure,  measured  aft  ot  apex 
XCP  center  of  pressure,  measured  aft  of  apex,  xcp 
Xcenl  area  centroid,  measured  aft  of  apex 

Subscripts 
LE  leading  edge 
TE  trailing  edge 

*  Member,  A1AA 
**  Associate  Fellow,  AIAA 


INTRODUCTION 

Typical  tactical  missiles  have  fins  with  small  chord 
lengths  relative  to  the  overall  body  length  Asa  resu  , 
static  stability  and  trim  are  not  significantly  affected  y 
errors  in  the  predicted  center  of  pressure  of  these  fins. 
Many  high-speed  vehicles,  however,  employ  very  highly 
swept  strakes  on  the  inboard  portion  of  the  wing.  In  some 
cases,  the  strakes  are  so  large  that  the  root  chord  of  the 
wing  is  as  long  as  the  fuselage.  The  stability  and  tom  of 
these  configurations  are  much  more  sensitive  to  errors 
the  predicted  wing  center  of  pressure. 

Preliminary  aerodynamic  design  codes  such  as  Missile 
Datcom  (Ref.  1)  and  the  Naval  Surface  Warfare  Center 
AP98  (Ref.  2)  use  approximate  methods  to  model  strake 
wings.  Given  a  stroked  wing  geometry,  they  define  an 
equivalent  straight  tapered  wing  which  is  used  for 
subsequent  analysis.  Spencer  (Ref.  3)  sh°wed  ^at  * 
equivalent  wing  method  can  produce  excellent  subsonic 
lift  predictions  if  the  appropriate  algorithm  is  used  to 
define  the  equivalent  straight  tapered  wing.  However,  he 
did  not  extend  his  analysis  to  include  supersonic  speeds  o 
consider  center  of  pressure  predictions. 

The  algorithm  used  for  multi-segment  wmg  centertf 
pressure  (XcP)  calculations  in  the  8/99  version  of  Missile 
Datcom  has  been  found  to  be  inadequate  whencornp^e 
with  data  from  future  generation  reusable  launch  vehicles 
with  long  strakes  such  as  the  X-34  and  X-37  s 
Figure  1  The  objective  of  this  effort  was  to  identify  an 
alternative  transformation  methodology  for  yieldi  g 
improved  xcp  predictions  for  multi-segment  wings  a 
small  angles  of  attack  suitable  for  inclusion  into  a  folure 
release  of  Missile  Datcom.  Only  planforms  with  ^break 
in  the  leading  edge  sweep  were  considered  P  anfo 
with  breaks  in  the  trailing  edge  will  be  studied 
subsequent  effort. 


American  Institute  of  Aeronautics  and  Astronautics 


A$C  -  0  2  ~  2  3  2  4 


JECHNICAI.  approach 

The  problem  of  developing  an  algorithm  for  multi- 

ST  rng  CI"ter  °f  pressure  calculations  can  be 
divided  into  two  d.stinct  parts:  (1)  defining  the  equivalent 

planfom  geometry  and  (2)  placing  the  center  of  pressure 
detennmed  from  the  equivalent  planform  onto  the  actual 
planform.  For  this  study,  nine  methods  for  defining  the 

mqeTn^nfP  anif0mi  ge°metry  Were  ana,yzed  and  three 
methods  for  placing  the  center  of  pressure  onto  the 

~Pla^  WT  COnsidered'  Detai,s  concerning 
each  of  the  methods  evaluated  as  part  of  this  investigation 

are  presented  next  followed  by  a  discussion  of  the 
approach  used  for  their  evaluation. 

Transformation  Methods 

Missile  Datcom  computes  the  center  of  pressure  of  a  wing 
using  design  charts  which  are  a  function  of  aspect  rat^ 

TTS  ?dmf  Cdge  SWCep  ang,e  and  Mach  ^mber.’ 
ah  of  the  transformation  methods  investigated  held 

(AR)  of  Ae  are.a  (?’/Pan  m)’  3nd  3Spect  rati0 
it.;,  A  °nglnaI  p,anform>  whiJc  allowing  other 
aspects  of  the  wing  geometry  to  vary.  For  purposes  of 

and  companson>  the  equations  used  for  each 
ot  the  transformation  methods  are  provided  in  Table  1  A 
more  formal  derivation  of  these  equations  follows.  ' 

Method  A 

The  wing  transformation  approach  used  for  center  of 
pressure  calculations  at  all  speeds  in  the  8/99  version  of 
Missile  Datcom  is  denoted  as  Method  A.  Missile  Datcom 

aL  m  *C.S  th,s  methodto  compute  the  normal  force  slope 
of  mult, -segment  wings  at  supersonic  speeds.  The 
method  evaluates  a  segmented  planform  by  transforming 
the  Planform  into  a  straight-tapered  approximation  with 

an  Prr ‘.at,°  “  **  °nginal  P|anform  and  defines 

weipht  r6  mgfedge  SWeep  anSle  based  on  the  area- 
weighted  average  of  the  leading  edge  sweep  angles  of 
each  segment,  using:  6 


Ale  (1) 

where  (Aif  j)  and  Sj  define  the  leading  edge  sweep  angle 
and  planform  area  of  each  wing  sejnenf,  respectively 
TTiejoot  chord  of  the  transformed  wing  is  calculated 


b(  \  +  X) 


Method  A  calculates  the  of  the  transformed  wing 
using  a  proportional  method  that  preserves  the  location  of 
the  xcp  as  a  percentage  of  the  root  chord  as  follows: 


ffl.  L[  xcp 


v  r  '  V  'effective 

Method  R 

Method  B  uses  the  transformation  developed  by  Spencer 
for  wings  at  subsonic  speeds.  Ref.  3.  Missile  So™ 

me?°d  t0  COmpUte  0,6  normal  force  slope  of 
ulti-segment  wings  at  subsonic  speeds.  It  compute!  an 

^  <*->  ■**  -  -a. 


Af/2=cos_l  -ZcosA„ 


This  is  then  used  to  calculate  the  effective  leading  edge 
sweep  angle  using:  5  ge 


A  tF  —  tan 


,  20-^) 
ar{\+X) 


A’  M*thod  B  the  taper  ratio  of  the 

g1  wing  and  uses  the  proportional  placement 
approach  provided  by  Equation  3  to  translate  the  X™  of 
the  transformed  fin  onto  the  original  fin.  ^ 

Method  C 

Method  C  is  a  variation  of  Method  B  that  uses  the  mid- 

V,r  I  7^7  a5  C  rather  ,han  cosine  of  the  angle  to 
calculate  the  effective  sweep: 


Kn  s^KnA  (6) 

Like  the  two  preceding  methods,  this  method  maintains 
the  taper  ratio  of  the  original  wing  and  locates  the  jl.  of 

EqU^^  P,anf0rm  proportiona,,y  wording  to 
Method  D 

iSJ*?  D  (RA  4)  differs  sig"ificant>y  from  the  previous 
methods  m  that  it  holds  constant  the  trailing  edge  sweep 

straiehMa  P  f***  i°f  ^  0rigina,  p,anform  t0  define  a 
aight-tapered  panel,  and  allows  the  taper  ratio  to  vary. 
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The  effective  leading  edge  sweep  angle  based  on  the 
trailing  edge  sweep,  aspect  ratio,  and  taper  ratio  of  t  e 
original  planform  is  calculated  using. 


A  ie  =  tan  1  tanAre+  — 


4  (1-4 

AR  (l  +  A) 


The  effective  root  chord  of  the  planform  is  computed 
using: 

,  (») 
Cr-  b  c, 

This  method  also  uses  a  different  procedure  to  translate 
the  Xcp  of  the  transformed  wing  back  onto  the  original 
planform.  Rather  than  placing  the  xcp  proportionally,  as 
the  previous  methods  do,  Method  D  lines  up  the  trailing 
edges  of  the  segmented  and  straight-tapered  planforms, 
and  places  the  xcp  of  the  original  planform  at  the  same 
distance  from  the  trailing  edge  as  the  xcp  of  the 
transformed  planform.  This  realignment  of  the  xcp  from 
the  transformed  wing  to  the  original  wing  can  be 
performed  using: 


i  =i-  i-p 

)  V  Cr  ) effective 


Cr,effective 


Method  E 

Method  E  is  a  variation  of  Method  D  that  maintains  the 
taper  ratio  of  the  original  planform  rather  than  the  tip 
chord.  It  also  uses  the  same  aft  XcP  trailing  edge 
alignment  method  as  Method  D  (Equation  9). 

Both  Methods  D  and  E  could  be  extended  to  planforms 
with  breaks  in  the  trailing  edge  sweep  angle  by  using: 


Ate  £  AjejS, 


Method  G 

The  equivalent  tangent  leading  edge  sweep  method, 
denoted  as  Method  G,  maintains  the  same  geometric 
properties  and  uses  the  same  Xq,  placement  method  as 
Method  A;  however,  it  uses  the  tangents  of  the  segmented 
wing  leading  edge  sweep  angles  to  determine  an  effective 
leading  edge  sweep  angle  using: 


A,p  =  tan  1  —  XtanAi^S,- 


Method  H 

Method  H,  obtained  from  Reference  2,  is  one  of  two  wing 
transformation  methods  contained  the  AP98  aerodynamic 
prediction  code.  The  method  holds  constant  both  the 
inboard  leading  edge  sweep  angle  and  the  trailing  edge 
sweep  angle  of  the  original  planform,  and  uses  the 
equations  below  to  calculate  new  root  and  tip  chords, 
respectively. 


rs  b, 

fs  b, 

nW 


Arp  ~  tan  A/jpj 


Are-tanAiEJ 


For  configurations  where  the  inboard  section  of  the  wing 
is  highly  swept,  such  as  with  strakes,  these  equations  can 
yield  negative  tip  chords.  This  method  was  not  evaluated 
where  a  negative  value  for  tip  chord  was  obtained  in  this 
study.  However,  AP98  assigns  a  tip  chord  of  zero  and 
calculates  a  new  root  chord  in  these  cases. 

The  Xcp  placement  approach  that  AP98  uses  lines  up  the 
area  centroids  of  the  original  and  transformed  wings  and 
places  the  XcP  of  the  segmented  wing  the  same  distance 
from  the  centroid  as  is  the  x^  of  the  straight-tapered 
wing,  using: 


Method  F 

A  hybrid  method,  denoted  as  Method  F,  combines 
approaches  used  in  Methods  A  and  D.  Like  Method  A, 
this  method  uses  Equation  1  to  calculate  an  effective 
leading  edge  sweep  angle  and  uses  the  proportional 
placement  approach  defined  in  Equation  2  for  locating 
Xcp.  However,  the  tip  chord  of  the  original  planform  is 
held  constant,  as  is  done  for  Method  D. 


'cp  |  _  (  xcenl  j  ,  Effective  x cp 


Method  I 

The  second  method  used  by  AP98,  denoted  as  Method  I, 
maintains  the  inboard  leading  edge  sweep  angle  of  the 
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original  planform  as  is  done  for  Method  H.  However, 
Method  I  maintains  the  taper  ratio  of  the  original 
planform  and  allows  the  trailing  edge  sweep  angle  to 
vary.  It  uses  the  area  centroid  placement  approach  used 
by  Method  H  to  locate  the  of  the  segmented  wing. 

Evaluation  Approach 

An  extensive  literature  search  was  performed  to  locate 
experimental  data  suitable  for  further  evaluation  and 
validation  of  the  most  promising  wing  transformation 
approaches.  Unfortunately,  very  little  data  for  isolated 
wings  is  available,  and  what  little  was  identified  exists 
only  at  subsonic  speeds.  In  order  to  perform  as 
systematic  evaluation  of  the  accuracy  of  the  various 
transformation  methods,  fourteen  segmented  planforms 
were  defined  to  represent  a  broad  range  of  reusable 
launch  vehicle  and  missile  fin  designs.  Referring  to 
Figure  2,  the  wing  geometric  parameters  varied  for  this 
study  included:  inboard  and  outboard  leading  edge  sweep 
angle,  trailing  edge  sweep  angle,  taper  ratio,  and  spanwise 
location  of  the  sweep  break.  The  span  of  the  planforms 
was  held  constant.  Table  2  lists  the  design  parameters  of 
the  segmented  planforms  investigated.  Corresponding 
scaled  drawings  of  these  wings  are  provided  in  Figure  3. 

Each  of  the  segmented  planform  designs  was  converted 
into  equivalent  straight  tapered  planforms  using  the  nine 
wing  transformation  methods  previously  described. 
These  configurations  were  analyzed  with  Missile  Datcom 
at  Mach  numbers  of  0.3, 2.0,  and  4.0  to  obtain  a  set  of  C,,, 
and  CN  predictions  from  which  an  Xcp  prediction  (defined 
in  terms  of  root  chord  percentage)  was  calculated  as: 

"7 L  =  ~~  (@  a  =  4  deg)  (15) 

Cr  ^ N 

In  order  to  deteimine  which  of  the  methods  examined  was 
most  accurate,  comparisons  were  made  between  Missile 
Datcom  results,  obtained  on  each  of  the  transformed 
straight  tapered  planforms,  and  “exact”  analytical 
solutions,  obtained  for  each  original  multi-segmented 
using  the  vortex  lattice  code  HASC  (Reference).  The 
HASC  analysis  was  performed  at  the  same  Mach  numbers 
and  angles  of  attack  as  the  Missile  Datcom  analysis.  The 
Missile  Datcom  results  for  all  of  the  transformed  wings 
assumed  proportional  placement  for  Xcp  (Eq.  3).  These 
were  post-processed  to  give  trailing  edge  placement  for 
Methods  D  and  E  and  centroid  placement  for  Methods  H 
and  I. 

To  address  numerical  uncertainty  issues  associated  with 
the  placement  and  density  of  the  bound  vortices  modeled 


in  HASC,  a  parametric  sensitivity  study  was  performed 
by  varying  the  number  of  chordwise  and  spanwise  bound 
vortices  for  a  limited  set  of  both  multi-sweep  and  straight- 
tapered  wing  geometries  for  which  experimental  data 
were  available.  Results  from  this  study  were  used  as  a 
guide  for  the  vortex  distributions  applied  on  each  of  the 
fourteen  multi-segmented  planforms  evaluated.  The 
distributions  used  for  assessing  the  Missile  Datcom 
results  were  sufficient  to  limit  numeric  difference 
associated  with  the  vortex  distributions  to  less  than  0.1% 
of  the  calculated  Xq/r^  value. 

As  an  additional  check  of  the  validity  of  using  HASC  to 
perform  this  evaluation,  three  of  the  multi-segment 
planforms  (Configurations  1,  3,  and  13)  were  analyzed 
using  the  space-marching  Euler  code  ZEUS  (Ref.  6). 
These  results  were  compared  with  the  HASC  results  for 
Xcp/Cf  at  Mach  2  and  4  and  found  to  be  within  1  percent  of 
one  another  at  a  four-degree  angle  of  attack. 

All.  predicted  results  obtained  for  this  investigation 
(Missile  Datcom,  HASC,  and  ZEUS)  were  independent  of 
Reynold’s  number  and  assumed  infinitely  thin  cross 
sections.  Sign  conventions  used  for  the  analysis  are 
shown  in  Figure  4. 

RESULTS  AND  ANALYSTS 

Table  3  provides  a  summary  of  the  average  error  in  center 
of  pressure  location,  as  a  fraction  of  the  root  chord  length, 
for  each  of  the  nine  methods  considered  at  Mach  0.3,  2.0 
and  4.0.  The  percent  error  is  defined  as  the  percent 
difference  between  the  results  obtained  using  the  various 
transformation  methods  and  the  “exact”  analytical 
solutions  obtained  from  HASC.  Methods  D  and  E  are 
clearly  superior  to  the  other  methods  with  average  errors 
on  the  order  of  V*  or  less  of  those  obtained  using  the 
transformation  method  used  in  the  8/99  version  of  Missile 
Datcom  (Method  A).  With  only  one  exception,  all  of  the 
methods  investigated  provided  greater  accuracy  than 
Method  A.  As  noted  above,  the  data  reduction  procedure 
generated  proportional  Xq,  placement  results  for  the 
planforms  generated  by  Methods  D,  E,  H  and  I.  These 
were  inferior  to  the  desired  placement  results  in  all  cases. 

Figure  5  plots  the  percent  error  in  wing  Xq,  versus  wing 
aspect  ratio  for  each  of  the  14  segmented  wings 
investigated  at  Mach  0.3,  2.0  and  4.0,  respectively.  In 
addition  to  the  results  obtained  from  Method  A,  results  for 
Methods  B,  C,  D  and  E  are  also  included.  The  number  of 
methods  shown  on  these  figures  was  limited  to  the  most 
promising  transformation  approaches,  in  order  to  facilitate 
clear  comparisons  between  them. 
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Inspection  of  Figure  4  shows  that  the  errors  associated 
lith  Method  A  are  frequently  greater  than  10  percen  at 
Mach  0.3  and  2.0  for  wing  aspect  ratios  less  than  .  • 
Methods  D  and  E  have  the  lowest  percent  errors  and  yield 
very  similar  results  for  all  configurations.  With  the 
exception  of  the  lowest  aspect  ratio  ^ngs  fconfi^rations 
10  and  1 1),  the  a*  percent  error  for  Methods  D .and  B  u 
less  than  5  percent  for  all  wings  at  both  Mach  0.3  and  2  . 

The  close  agreement  in  results  is  not  surprising  given  the 
similarity  of  the  two  methods.  At  Mach  4.0,  all  metho 
Z2a,ed  gre,w  accttracy,  but  the  ever*! bentfe 
observed  remained  consistent  with  those  observed  at 
Mach  0.3  and  2.0.  The  xcp  percent  error  for  Methods  D 
and  E  were  both  less  than  5  percent  for  all  configurations 
evaluated  at  Mach  4.0. 

Methods  A  and  E  were  compared  with  a  limited  set  of 
subsonic  experimental  data  obtained  References  7  through 
H  to  validate  the  assessments  made  from  cot^sot* 
with  “exact”  solutions  generated  by  HASC.  Tab 
identifies  the  wing  configurations  and  test  condemns  fo 
which  subsonic  experimental  data  was  obtained  for 
validation  purposes.  Figure  6  provides >  the scakd 
planform  views  of  each  of  the  wmgs  identified  m  Table  4. 

A  plot  comparing  xcp  percent  error  as  a  function  of  aspec 
ratio  for  both  Methods  A  and  E  is  shown  m  Figure? .A. 
observed  in  “exact”  data  compansons  made  m  Figu  , 
Method  E  XcP  predictions  are  consistently  wi  in 
percent  of  the  experimental  data,  while  Method  A  Xcp 
errors  are  anywhere  from  two  to  five  times  greater  than 
Method  E  The  maximum  error  for  the  expenmenta  Xcp 
data  is  estimated  to  be  within  one  percent  of  the  wing  root 
chord. 

The  values  of  center  of  pressure  discussed  thus  far  are 
only  used  by  Missile  Datcom  for  the  linear  lift 
contribution,  i.e„  at  low  angles  of  attack^  Miss* i  Dateom 
uses  the  wing  centroid  as  the  center  of  presme  for  the 
non-linear  lift  contribution.  The  overall  wing 
computed  using; 

xcp  CNJinear  (Xcp  lcr)  +  CN, non- linear  (* cen I  [  Cr)  (J6) 

Cr  ~  CN 

To  assess  the  performance  of  Method  E  relative  to 
Method  A  in  Missile  Datcom,  plots  of  Xq/Cr  were 
generated  versus  angle  of  attack  for  wing  configurations 
IS  16  17  22  and  23.  These  plots  are  presented  m 
Figures  8  through  12  in  numeric  order.  Method  E^eids 
significantly  better  yields  results  than  Method  A  for  a 
four  configurations  and  Method  E  is  within  four  percent 


of  experimental  x^/c,  values  for  angles  of  attack  up  to  25 
degrees. 

grnvfMARY  AND  CONCLUSIONS 

The  results  of  this  study  have  shown  that  several 
transformation  methods  are  superior  to  the  method 
contained  in  the  8/99  version  of  Missile  Dateoire  The 
independent  validation  of  findings  from  the  exact 
solutions  study,  obtained  through  the  use  °f  experanental 
data  has  confirmed  that  the  incorporation  of  Method  E 
into  Missile  Datcom  will  result  in  significant 
improvements  to  the  capability  of  the  code  to  predict  the 
-  of  multi-segmented  wings  at  subsonic  speeds.  It  is 
also  reasonable  to  assume  that,  based  on  the  results  of  the 
“exact”  study.  Method  E  will  also  yield  significant 
improvements  for  supersonic  predictions. 

Due  to  the  inability  to  predict  “exact”  analytical  solutions 
and  the  lack  of  experimental  data,  no  conclusions  can  be 
drawn  for  the  transonic  flight  regime.  There  is  clearly  a 
need  to  obtain  experimental  data  for  validating  e 
at  supersonic  transonic  speeds. 

The  next  release  of  Missile  Datcom  (Version  9/02)  will 
incorporate  Method  E  for  predicting  the  x^  of  multi- 
segmented  wings.  This  improvement  will  be  of  particular 
benefit  for  the  aerodynamic  assessment  of  reusable  launc 
vehicle  design  concepts  as  well  as  unconventional  missile 
fin  designs. 
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TaWe  K  Summary  of  Transformation  Methods  for  Multi-Segment  Wings 
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Table  2.  Fin  Geometries  for  Systematic  Investigation 


Fin 

A  le, 

i 

A  LE.2 

Ate 

■ 

R 

AR 

1 

80 

0 

0 

0.25 

0.5 

1.22 

2 

80 

40 

-40 

0.25 

0.5 

0.72 

3 

80 

40 

-20 

0.25 

0.5 

0.84 

a 

80 

40 

0 
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0.5 

0.94 

80 

40 

20 

0.25 

0.5 

1.10 

6 

80 

40 

40 

0.25 

0.5 

1.42 

a 

80 

40 
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1.28 

8 

80 

40 
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40 

0 

0 
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3.82 

13 

0 

40 

0 

0.5 

0.5 

2.72 

14 

0 

40 

0 

0.25 

0.25 

3.32 

Table  3.  Average  Error  in  Predicted  xcp 


Method 

M  =  0.3 

M  =  2.0 

M  =  4.0 

A 

12.2 

11.0 

7.5 

B 

5.9 

6.0 

3.3 

C 

8.6 

7.8 

4.2 

D 

2.6 

3.3 

1.8 

E 

2.3 

2.8 

1.1 

F 

8.2 

8.1 

6,0 

G 

13.5 

12.3 

13.3 

H 

9.5 

9.6 

6.0 

I 

7.7 

9.3 

6.0 

Table  4.  Fin  Geometries  Used  for  Validation 
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2 
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9 

22 

76 
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0 
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0.23 
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23 

80 
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Figure  I.  Photographs  of  X-34  and  X-37. 


Figure  4.  Sign  Conventions. 


Figure  2.  Wing  Parameters. 
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Figure  5.  Analytical  Error  Versus  Aspect  Ratio. 
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XCP  (PERCENT  ROOT  CHORD) 
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Figure  6.  Wing  Planforms  Used  for  Validation. 


ASPECT  RATIO 


Figure  7.  Experimental  Error  Versus  Aspect  Ratio. 


ANGLE  OF  ATTACK  (DEG) 

Figure  10.  Wing  Xq,  Versus  Angle  of  Attack, 
Configuration  17,  Mach  =  0.14. 
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Figure  8.  Wing  xq,  Versus  Angle  of  Attack, 
Configuration  15,  Mach  =  0.14. 
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Figure  12.  Wing  xcp  Versus  Angleof  Attack, 
Configuration  22,  Mach  =  0.23. 
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